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In this paper we investigate isospin transport mechanisms in semi-peripheral collisions at Fermi 
energies. The effects of the formation of a low density region (neck) between the two reaction partners 
and of pre-equilibrium emission on the dynamics of isospin equilibration are carefully analyzed. We 
clearly identify two main contributions to the isospin transport: isospin diffusion due to the N/Z 
ratio and isospin drift due to the density gradients. Both effects are sensitive to the symmetry 
part of the nuclear Equation of State (EOS), in particular to the value and slope around saturation 
density. 
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I. INTRODUCTION 

In the last few years the increased accuracy of the ex- 
perimental techniques has renewed interest in nuclear 
reactions at Fermi energies. Exclusive measurements, 
event-by-event analysis, and a 4-7T coverage allow a deeper 
investigation of the evolution of the reaction mechanisms 
with beam energy and centrality. New insights into the 
understanding of the nuclear matter equation-of-state 
(EOS) were gained pj . In particular, recent experimen- 
tal and theoretical analyses were devoted to the study of 
the properties and effects of the symmetry term of the 
EOS (asy-EOS) away from saturation conditions 

Indeed, the two-component character of nuclear mat- 
ter adds some special interest to the dynamics of heavy 
ion collisions at intermediate energies, between 20 and 
lOOAMeV. In central collisions isospin distillation is an 
important effect in multifragmentation of charge asym- 
metric systems. Here phase separation is driven by 
isoscalar-like unstable fluctuations, i.e. local in phase 
variations of proton and neutron densities 0, IE S 0- 
This leads to a more symmetric "liquid" phase of frag- 
ments surrounded by a more neutron rich "gas" relative 
to the original asymmetry of the system. Isoscaling phe- 
nomena, observed experimentally, provide indications for 
such a scenario H, @ . 

In semi- peripheral collisions between nuclei with dif- 
ferent N/Z ratio, isospin dynamics will drive the sys- 
tem toward a uniform asymmetry distribution. The de- 
gree of equilibration, correlated to the interaction time, 
should provide some i nsig hts into transport properties of 
fermionic systems [l(j, [Tlj , in particular give information 
on transport coefficients of asymmetric nuclear matter 

mm. 

The aim of this work is to investigate the isospin trans- 
fer through the neck region in semi-peripheral collisions 



of asymmetric nuclei at Fermi energies. Of particular 
interest is the role of the density dependence of the sym- 
metry energy in this process. The isospin transfer was 
measured for collisions of different Sn isotopes at MSU 
[Til flEl fl6| and interpreted theoretically with the result 
that the a.sy-EOS should be rather stiff. In these works 
the effect of pre-equilibrium emission, which changes the 
isospin content of the interacting system, was not ana- 
lyzed explicitly. Here we will discuss this questions in 
detail, as well as the different transport processes affect- 
ing the final isospin content. Quantitatively, dynamical 
isospin effects can be properly understood only from mi- 
croscopic calculations based on transport models. We 
will base our study on a stochastic BNV transport model 
(seerefs. 0113 for more details on the main ingredients 
of this approach). 



II. ISOSPIN EQUILIBRATION PROCESS 

We are focusing on the charge asymmetric collision 
124 Sn+ 112 Sn, at bOAMeV bombarding energy, to which 
we refer as the mixed system, (M). To investigate 
the density (p) dependence we consider here two rep- 
resentative parameterizations of the symmetry energy, 
E sym (p,I)/A = C sym (p)I 2 ,I = (N - Z)/A : one show- 
ing a rapidly increasing behaviour with density, roughly 
proportional to p 2 (asysuperstiff) and one where a kind 
of saturation is observed above normal density (soft, 
SKM*) (see Ref.[l|| for more detail). 

The BNV simulations have been performed for 
semi-peripheral collisions at impact parameters b = 
6,8,9, 10/m. In the last two cases the reaction has 
dominantly a binary character and the charge asymme- 
try of primary projectile (target)- like fragments, PLF 
( TLF), should provide the essential information about 
the isospin equilibration rate. At b = 8/m already about 
25% of the events are ternary. An IMF can be formed 
in the mid- velocity region by neck fragmentation |l9j| . 
For more central events this mechanism becomes dom- 
inant: at b = 6fm, one or two IMF's are found in 
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FIG. 1: 124 Sn + 112 Sn b = 8/m (upper 7muj and 6 = 10/m 
(lower row) collision: density contour plots. The side of each 
box is 40 fm. 

more than 70% of events @. The fragment formation 
in the neck region will influence the final isospin distri- 
bution of the PLF/TLF, and will render considerably 
more difficult the interpretation of the results. Thus in 
the following we select only binary events in our anal- 
ysis. We define the average interaction time, t c , as the 
time elapsed between the initial touching and the mo- 
ment when PLF and TLF reseparate. From our simula- 
tions we obtain t c m 140, 120, 100, 80 fm/c for the impact 
parameters b = 6,8,9, 10 fm, respectively. Four hundred 
events were calculated for each initial condition and for 
each asy-EOS. 

Typical density contour plots, at b = 8/m and b — 
10/m, are shown in Figure ^ We note the dynami- 
cal evolution of the overlap region: driven by the fast 
leading motion of PL— and TL— prefragments, the for- 
mation of a lower density interface can be clearly ob- 
served after around 40/m/c. An isospin migration, or 
transport, takes place during this transient configuration 
of two residues with densities close to the normal one, 
separated by a dilute neck region. In contrast, in deep- 
inelastic collisions at lower energies the isospin equilibra- 
tion is driven by the N/Z difference between the interact- 
ing nuclei having a quite uniform density profile without 
a low density interface until separation [2fJ . 

We quantify the degree of equilibration by the isospin 
transport (imbalance) ratio [2l|, defined as: 
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Here i = P,T stands for the projectile-like (target- 
like) fragment. The quantities Ii refer to the isospin, 
/ = [N — Z)/A, or in general to any isospin depen- 
dent quantity, characterizing the fragments at separation 
time, for the mixed reaction (M, 124+112), the reactions 
between neutron rich (H, 124 + 124), and between neu- 
tron poor nuclei (L, 112 + 112), respectively. A value of 
Ri approaching zero is an indication of a larger degree of 
equilibration. The extreme cases Rp = — 1 and Rp = 1 
correspond to the absence of any isospin transfer. 
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FIG. 2: The isospin transport ratio Rp,t, eg. (1), for the 
asysoft (full squares) and asysuperstiff (circles) EOS's as a 
function of the interaction time t c , corresponding to different 
impact parameter b. The band betwee n the tw o solid lines 
corresponds to the experimental data of JflA. 



The results of the calculations are shown in Figure [21 
for the dependence of Rp/p on the interaction time t c for 
the asysoft (squares) and asysuperstiff (circles) EOS's. 
The figure also shows the experimental values extracted 
in ref . [TEl IT^ at semi-peripheral collisions at about b = 
8/m. We conclude, as in ref. [l^, that an asystiff-likc 
EOS provides a better agreement with the experimental 
observations. 

A significant difference between the two equations of 
state is evident for larger interaction times, i.e. smaller 
values of the impact parameter, 6 = 6, 8/m. The smaller 
values of isospin transport ratios for the asysoft EOS 
point toward a faster equilibration rate. In refs. 0, |2^ 
an explanation was based on the observation that below 
normal density the asysoft EOS has a larger value of the 
symmetry energy. Therefore an enhanced isospin equili- 
bration will occur if the diffusion takes place at uniform 
lower density. We intend to show, that in fact the mech- 
anism of charge equilibration is more complicated due 
to dynamical evolution of the reaction at these energies. 
Fast particle emission and density gradients will also play 
a role. In the next section we investigate more in detail 
the various influences on the isospin transfer process. 



III. ISOSPIN SHARING AT FERMI ENERGIES 

The isospin content of the two residues in a mixed col- 
lision system at separation time is determined by the 
interplay between the particle emission to the gas from 
each nucleus during the overlap and the transfer of nu- 
cleons through the neck. We thus write simple balance 
equations: 
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Here I P (I T ) and A P (At) are the PLF (TLF) asym- 
metry and mass at separation, I p (Ij,) and A° p (A 1 ^) the 
initial projectile (target) asymmetry and mass. Then 
IgP (Igr), A g p (AgT) are the asymmetries and masses 
of the projectile/target "gas", i.e. of the pre-equilibrium 
particles emitted by the projectile (target) during the in- 
teraction time. Finally Ipt (Itp), and Apt (Atp) are 
the asymmetry and mass of all nucleons transferred from 
projectile (target) to target (projectile). 

In Figure|3]we plot the time evolution of the quantities 
IgP , IgT and A g p , A g T , as well as the values of Ipt and 
Itp for the asysoft and the asysuperstiff EOS for two im- 
pact parameters. We note that I g p is much larger than 
I P . The same is true for the target but the difference is 
smaller. Thus the pre-equilibrium emission reduces the 
N/Z difference between the two nuclei, competing with 
the transfer process. This is also clearly seen in Figure 
0] where the time evolution of the projectile, target, and 
composite system apparent asymmetry, I p pp , I T PP and 
7^? p is shown. The apparent asymmetry is calculated 
taking into account only the preequilibrium emission but 
not the isospin transfer inside the matter. A stronger 
variation is observed for the neutron rich projectile. The 
composite system asymmetry, I^ pp , is the limiting value 
toward which the two participants would evolve in com- 
plete isospin equilibrium. In contrast, we show in the fig- 
ure the final asymmetries Ip and It including the trans- 
fer processes. 

Comparing the results for the two EOS's in Figs. 3 
and 4 we see that a more neutron rich composition of 
pre-equilibrium emission is generated in the asysoft case 
because below normal density, from where most of the 
emitted nucleons originate, the neutrons (protons) are 
less (more) bound than for the asysuperstiff EOS. The 
differences between the two asy--EOS"s are strongly re- 
duced at larger impact parameters, as seen in the results 
for b = 10 fm in Figures |3] and 0] since the interaction 
times are much shorter. 

We next show in Figure |S] the dependence of the quan- 
tities introduced above on interaction time, i.e. impact 
parameter. Here the projectile (target) and gas final 
asymmetries Ip ,I g p (It, I g r) are confronted with the 
transferred asymmetries Ipt and Itp- For the projectile, 
that is neutron rich, both pre-equilibrium emission and 
nucleon transfer drive the system toward a more sym- 
metric configuration (left). The two processes have op- 
posite effects and thus tend to compensate for the tar- 
get (right). Therefore the projectile asymmetry has a 
more pronounced deviation from the corresponding ini- 
tial value in comparison to the target. 

We find a a clear dependence on the asy-EOS of the 
isospin transferred between the two nuclei. In particular, 
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FIG. 3: 124 Sn + 112 Sn collit 
b = 10/m (right): Time evolution of isospin content and 
mass of preequilibrium particles l g p, A g p (I 9 t, A 9 t) emit- 
ted from the projectile (top panels) (resp. target, lower pan- 
els) for asysoft (solid lines) and asysuperstiff (dashed lines) 
EOS. The squares (circles) indicate the mass and isospin of 
the nucleons that have migrated from projectile to target Ipt 
(and vice versa, Itp ) for asysoft (full symbols) and asysuper- 
stiff (empty symbols) at separation time. The horizontal lines 
in the left panels give the initial isospin of the projectile, Ip, 
and target,I%, respectively. 
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FIG. 4: 124 Sn + 112 Sn at b = 8/m (left windows) and 
b = 10/m (right windows): Influence of pre-equilibrium emis- 
sion on isospin content of projectile (top) and target (bot- 
tom) systems. The curves show the "apparent" isospin with- 
out transfer, while the symbols denote the real isospin at sep- 
aration time including also tranfer. Dashed lines and open 
symbols correspond to asysuperstiff EOS and the solid lines 
and full symbols to the asysoft case. 
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FIG. 5: 124 Sn+ 112 Sn collision with asysoft (top) and asysu- 
perstiff (bottom) EOS as a function of collision time or im- 
pact parameter. Left panels: Ip (black circles), I g p (crosses), 
Ipt (rhombs) and Itp (squares). Right windows: It (black 
circles), I b t (crosses), Ipt (rhombs) and Itp (squares). The 
horizontal solid lines refer to the initial asymmetry of the pro- 
jectile (target). 



for b = 8fm(t c = 120/m/c) we observe that 

j-(asysuperstiff) ^ j(asysoft) > j0 _ q 
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We will show in the next chapter that the origin of 
these inequalities lies in the existence of the low den- 
sity interface and the density dependence of symmetry 
energy. The asysuperstiff EOS favors the neutron mi- 
gration toward the neck region from both participants. 
This explains why simultaneously j^y su P erstl f i) > jo 



and /, 



(asy super stif /) 



TP 



PT ^ ± P 

> I T . We will see that for the asysoft 



EOS this effect is weakened. 



IV. INTERPRETATION OF THE RESULTS 

The above arguments can be made more explicit by 
considering that the proton and neutron migration is dic- 
tated by the spatial gradients of the corresponding chem- 
ical potentials p p / n {p p , p n ,T) j2^|. The currents of the 
two species can be expressed as follows 



jn = -CtVfl n (p p ,pn,T) 



dp r , 
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where ct is a constant. Rewriting these expressions in 
terms of Landau parameters, (|j| and references 

therein), using: 
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(where N q (T) is the level density), we obtain: 
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U q is the neutron or proton mean-field potential. The 
second lines in eqs. (7,8) define p and n drift and diffusion 
coefficients due to density and isospin gradients, D? and 
Dg, respectively. The terms N~ 1 (1±I) and fP ± lf q , 
that appear in the density drift coefficients DP, can be 
expressed as 



N~ 1 (l±l) = 2N- 1 {p,T)±4I Jj ym 



0(J 2 ); 



Qnpot 

fP q ±lf q =F{p)±AI-^ 



+ 0(I 2 ),(+n,-p), (11) 



where the function F(p) depends only on the isoscalar 
part of the interaction. We have denoted by and 
C v s ° l m the kinetic and the potential part of the symme- 
try energy coefficient C sym , respectively. Combining eqs. 
(|llf> . one can see that the isovector part of the nuclear 
interaction enters the coefficients DP through the deriva- 
tive of the total symmetry energy C sym . On the other 
hand the isospin diffusion coefficients D 1 are proportional 
to the quantity 



PiNq 1 + fq) 



A[C sym ±I(p^^ 



Csym)], (12) 



and depend, in leading order, on the value of the sym- 
metry energy coeffcient C sym . 

Various particular situations can be derived from these 
relations. In symmetric nuclear matter D* n = —D p and 
DP = D^. In the absence of density gradients the proton 
current will flow oppositely and equal in magnitude to 
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FIG. 6: Ratios of drift coefficients R q 

r^i.asy super sti f f / j-)i,asy so ft ■ „ j 7 „ „ „ „„ 

D q " y JJ jU q " J ,i = p, 1 and q = n,p, as 
function of the density for fixed asymmetry I — 0.2. 



as many neutrons are removed as for the asysoft EOS. 
Also, as it was shown, for the asysuperstiff EOS more 
asymmetric matter is also transferred from the target. 

We now estimate the effect of isospin transport and 
pre-equilibrium emission on the transport ratios. Within 
the following approximations 
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where the indices M, 77, L refer again to the mixed 
(124+112), the neutron rich (124 + 124) and the neu- 
tron deficient (112 + 112) systems, we arrive at a sim- 
plified expression for the isospin transport ratio for the 
projectile which shows explicitely the dependence on the 
isospin transport (Ipt — Itp) and on the pre-equilibrium 



emission (I^ P 



Igp)-- 



the neutron current. On the other hand, for density gra- 
dients only, in asymmetric nuclear matter the proton and 
neutron currents may have the same direction but assume 
different values, inducing isospin gradients (24|. Such a 
situation can be encountered in semi-pheripheral colli- 
sions between identical, charge asymmetric nuclei with 
the formation of a dilute intermediate region. 

For the two asy-7?OS"s we calculate the coefficients 



p,I and q — n,p. We plot the ratios Rl = 



D i,asy super stiff j D i,asysoft in Figure | as a function of 
the density for a fixed asymmetry 7 = 0.2 . These values 
of the asymmetry and density are close to the physical 
conditions expected for the projectile or target region. 
The only negative coefficient is D^. The isospin gradi- 
ents, directed from the projectile to the neck and from 
the neck to the target, induce neutron and proton flows 
in opposite directions. However the ratios of the corre- 
sponding coefficients are quite close to unity for the two 
asy-EOS's and therefore the effects are similar. 

Since the density gradient is oriented from pro- 
jectile and target residues to the neck, neutrons 
and protons migrate from higher toward lower den- 
sity regions. Around and below saturation density 

RP = D P,asysupersUff / D p,asysoft > j and R p = 

D p,asysu P erstiff/ D p,asysoft < 1 These inequalities sug- 
gest that more neutrons and less protons migrate from 
projectile toward neck in the case of asysuperstiff EOS 
resulting in the formation of a more neutron rich inter- 
mediate region. This is due to the larger value of the 
derivative of the symmetry energy around normal den- 
sity and is in agreement with the behavior observed in 
the numerical simulations, see eqs. (4|SJ. 

We note that for asysuperstiff EOS, in spite of an en- 
hanced isospin migration toward the neck at separation 
time, the projectile residue is more asymmetric in com- 
parison to the asysoft case. One reason is that for the asy- 
superstiff EOS during the pre-equilibrium emission not 
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With similar approximations the target isospin trans- 
port ratio can be expressed as: 
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It is observed that the transport ratios depend on the 
difference Ipt ~ Itp as expected. However, in contrast 
to what was assumed in refs. 0, 12^, it is seen that 
they also depend on the pre-equilibrium emission which 
reduces the absolute value of the transport ratios. Both 
effects are smaller in the case of the asysuperstiff EOS, 
for which indeed a larger R ratio is obtained. 



V. CONCLUSIONS 

In this work we have studied processes related to 
isospin equilibration in semipheripheral collisions at 
Fermi energies and their dependence on the symmetry 
term of the EOS. A special feature of these reactions is 
the development of a low density interface between the 
two residues. The neck region is controlling the proton 
and neutron currents and their direction. The presence 
of density gradients also affects the isospin exchange be- 
tween projectile and target and we have shown that this 
is sensitive to the density dependence of the symmetry 
energy. The neutron to proton ratio emitted during the 
interaction stage is also influenced by the asy- EOS. The 
interplay between the two processes leads to a stronger 
equilibration for asy-soft EOS, as it is evidenced by the 
isospin transport (imbalance) ratio. Actually, in the asy- 
stiff larger isospin transfer is observed, due to 



6 



the presence of density gradients, directed from PLF 
and TLF towards the neck region. However, since we 
are studying binary processes, finally we observe a kind 
of compensation between the asymmetry of the matter 
transferred from projectile to target (Ipr) and from tar- 
get to projectile(Trp)- /.From this point of view, to put 
in better evidence effects due to the presence of density 
gradients, it would be more appropriate to study events 
where fragments originating from the neck region are also 
detected, possibly with their isospin content. 

In the present study a rapidly increasing symmetry en- 
ergy at subnormal densities (asysuperstiff) appears to be 
in better agreement with the existing data [lq. More 
recent work also considered momentum dependent inter- 



actions, showing that the more repulsive character of the 
overall dynamics may reduce the symmetry energy stiff- 
ness required to reproduce the data 22] . The study of 
the interplay between the effects due to the isoscalar and 
isovector part of the EOS on isospin transport observ- 
ablcs deserves further attention. 

In conclusion, charge equilibration measurements in 
semi-peripheral heavy ion collisions at Fermi energies 
provide new independent observables to study the poorly 
known density dependence of the symmetry term of the 
nuclear EOS. This is of interest for other properties of 
asymmetric matter, like neutron skin and isovector col- 
lective response in finite nuclei, and may also be impor- 
tant for neutron star crust structures yl. 
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